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ABSTRACT 

The excursion set theory of halo formation is modified by adopting the fractional Brow- 
nian motion, to account for possible correlation between merging steps. We worked 
out analytically the conditional mass function, halo merging rate and formation time 
distribution in the spherical collapse model. We also developed an approximation for 
the ellipsoidal collapse model and applied it to the calculation of the conditional mass 
function and the halo formation time distribution. For models in which the steps are 
positively correlated, the halo merger rate is enhanced when the accreted mass is less 
than ~ 25M*, while for the negatively correlated case this rate is reduced. Compared 
with the standard model in which the steps are uncorrelated, the models with posi- 
tively correlated steps produce more aged population in small mass halos and more 
younger population in large mass halos, while for the models with negatively correlated 
steps the opposite is true. An examination of simulation results shows that a weakly 
positive correlation between successive merging steps appears to fit best. We have 
also found a systematic effect in the measured mass function due to the finite volume 
of simulations. In future work, this will be included in the halo model to accurately 
predict the three point correlation function estimated from simulations. 

Key words: cosmology: theory - large scale structure of the Universe - galaxies : 
halos - methods : analytical 



1 INTRODUCTION 



then satisfies a diffusion equation 



The excursion set theory provides a simple and intuitive 
model for cosmic structure formation. In this theory, as one 
varies the smoothing scale R, the linear density fluctuation 
8 a obtained with smoothing window W(R) forms a one di- 
mensional random walk. It is deemed that the non-linear 
halo formation and evolution history can be treated with the 
corresponding excursion set theory: a halo is formed when 
a pre-set barrier S c (z) is exceeded by a tracjectory of the 
random walk. 

It is convinient to take the variance of the random den- 
sity field S(R) — cr 2 (R) as the pseudo-time variable. Prop- 
erties of the random walk rely on the window function used 
and the nature of the primordial fluctuation. For Gaussian 
fluctuations and using the sharp k-space filter, the random 
walk produced by smoothing is a normal Brownian motion 
(hereafter NBM), i.e. there is no correlation between steps. 
The density of trajectories Q(S,5) passing through (S,8) 
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From Q(S,S) one could derive the halo mass function a nd 
merger rates, etc. (jBond et al.lll99ll : lLacev fc Cordll993l ). 

However, although the N-body simulation results gen- 
erally agree with the predictions of the excursion set the- 
ory, there are significant deviations in the details. The dis- 
crepancy is particularly severe in the description of small 
mass halos. This discrepancy can be partly overcomed 
by replacing the spherical collapse model with ellisoidal 
collapse model, i.e. by introducing a moving barrier in- 
stead of a fixed barrier imposed on the random walk (e.g. 
ISheth fc Tormenl l2002h . This practical approach provides 
a reasonably good fitting formula for halo mass function 
l|Warren et al. 20061 ). but its prediction on the formation 
time distributi on of low- and intem ediate-mass halo remains 
unsatisfactory jGiocoli et al.|[2007l) . One suspects that if cor- 
relations between steps of the random walk is introduced, 
the excursion set theory might be improved. Indeed, if the 
smoothing window function is not a sharp k-space filter but 
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a Gaussian filter or a real space tophat filter, the excursion 
steps would be correlated. 

In a previous work by the leading author, the fractional 
Brownian motion (FBM), the simplest random walk with 
steps correlated in long range, was introduced to generalize 
the excursion set theory. The correlation between steps of 
the random walk was shown to be capable of m odifying th e 
final halo mass function in a non-trivial way 
The model presented there was incomplete though, as the 
solution for ellipsoidal collapses and the treatment of halo 
merging history were not discussed. These are the topics of 
the present report. 

This paper is organized as follows. In section 2 we fit the 
FBM into the excursion set theory to account for the possi- 
ble correlated halo formation process. Then in section 3 the 
diffusion equation of FBM is solved in the spherical collapse 
model, and we calculate the conditional mass function, halo 
merger rate and the halo formation time distribution with 
the modified theory. Treatment to the ellipsoidal collapse 
is given in section 4, together with a comparison between 
theoretical predictions and measurements with simulations. 
The final section contains our conclusions and discussion. 
We adopt the flat ACDM model with the following set of 
cosmological parameter values: Q m — 0.3, Qa = 0.7, h — 0.7 
and as = 0.9. 



2 MODELING CORRELATED MERGING 
STEPS 

2.1 The diffusion equation 

To work out the halo conditional mass function, the central 
element is the diffusion equation which governs the behavior 
of the random walk. Such diffusion equation in turn depends 
on the understanding of the random walk with which the 
physical problem is concerned with. The conditional mass 
function is in fact a two-barriers crossing problem of a ran- 
dom walk, which means the key object we shall check is the 
scaling relation between 8(Si) — S(So) and Si — So. 

The FBM is a random process X(t) on some probability 
space such that: 

(i) X(t) is continuous and X(0) = 0; 

(ii) for any t > and r > 0, the increment X(t + r) — X(t) 
follows a normal distribution with mean zero and variance 
r 2a , so that 



P(X(t + r)-X(t) <x) = 



2 /2t 



du . (2) 



The parameter a is the Hurst exponent, if a = 1/2, it i s 
reduced to the normal Brownian motion (c.f. iFederl [l988h . 
It is also easy to see that the following is satisfied by the 
FBM: 



[X(t- 



X(t)} 2 ) 



{X(t)[X(t + r)-X(t)]) = 



(t + r) 2a - t 2a - r 2a 



(3) 



With this definition of FBM, the trajectory density Q a (X, t) 
at time t in interval (X, X + dX) follows the diffusion equa- 
tion fc.f. lLutd[200lh . 

^r^ v ^x 2 -- v = 2d- t {x{t)) = at ■ (4) 



Why invoke FBM? The trajectory of S(S) is character- 
ized by properties of the increments S(Si) — S(So) between 
any pairs of (So, Si > So)- If the smoothing window func- 
tion corresponding to halo definition is not a sharp k-space 
filter but e.g. a Gaussian or top-hat, it is easy to check 
that ([5(Si) — S(So)] 2 } is not Si — So but proportional to 
(Si — So) Q with a constant a in broad scale range. The scal- 
ing relation is the flagging attribute of the FBlvQ, so it is 
reasonable to install the FBM into the excursion set theory 
as the simplest approach to analytically inspect halo models 
constructed with random walks with correlated steps. 

On the other hand, we know in practice that the bound- 
ary, and subsequently the mass of a halo identified in a 
simulation is rather arbitrary, e.g. the mass picked up is 
often defined by a halo's virial radius while the mass out- 
side is simply not counted. The S from the virial mass M 
is very likely not the real place where the particular ran- 
dom walk correspon ding to the halo hit the barrier. In fact 
ICuesta et al.l (|2007t ) argued that if one replaces the virial 
mass with the static mass, the mass function agrees with 
the Press-Schechter formula remarkably well, rather than 
the Sheth-Tormen one, and the ratio of the virial mass to 
the static mass depends on redshift and halo mass. 

While the true scale can be smaller or larger than the 
S inferred from R(M), one can always parametrize the true 
scale with S(R) so that {8 2 } = S 2a . Taking the approxima- 
tion that a is a constant not too far away from 1/2, we can 
comfortably assume the applicability of the FBM. 

For an FBM with (S 2 ) = S 2a , by Eq. (© , the variance 
of the increments at two points (So, Si > So > 0) satisfies 



([<5(S 1 )-<5(So)] 2 )=(Si-So) 2 



with a € (0, 1) , Si > So > 



(5) 



so that a new trajectory is formed by 5(AS) = S(Si) — 
5(So) along AS = Si — So, given any source point (So,^o = 
5o(So)). Apparently the new random walk also complies with 
the definition of FBM, its diffusion equation being 



dQ a _ d 2 Q 
OAS ds 2 



T>AS 



1 d(S 2 ) 

2 dAS 
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It is this diffusion equation that ought to be solved to fig- 
ure out the conditional mass function, which can be fur- 
ther transformed to the familiar diffusion equation of normal 
Brownian motion 

dQ a 1 d 2 Q a 



os 

by the substitution S = AS 2 



2 df 



(7) 



2.2 The correlation 

Imagine a halo is formed by the collapsing condition <5(Si) = 
lui at Si = a 2 (Aii) with Mi being the halo mass. After 
some time the halo merges into a bigger halo with mass 
Mo > Mi at So = cr 2 (Mo) < Si by another collapsing con- 
dition <5(So) = ujq < Wi. The following question needs to be 



There are many anoma lous random walks such as the fractal 
time process jLutd l200ll) which have the same scaling feature 
and are classified as sub-diffusion. The FBM is the simplest one 
of them. 
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Figure 1. Correlation functions between successive halo merging 
steps (Eq|9]l. 



addressed: how the formation event by uio at So is correlated 
with the past merging process of ivi — loq within Si — So? 

If this process is approximated by the FBM, the cor- 
relation function can be easily calculated with Eq. <j3j . We 
have 

£0,0^1 = ((wo - 0)(wi - uio)) 



s 2a 

2 



£0 
Si 



(8) 



It is more instructive to capture the correlation between 
two successive merging steps: U12 — wi within S2 — Si, and 
oji — u)q within Si — So (S2 > Si > So, 1^2 > tui > u>o). The 
correlation function is similarly 

(S 2 - S ) 2a 



£0- 



»l,l->2 = ((^2 — Ct>l)(wi — Wo)) = 

Si — So 



x 
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which is positive when a > 1/2 and negative if a < 1/2. 
Apparently, a — 1/2 causes null correlation (Fig. [1]). Thus 
a unified paradigm is given by a single parameter controlled 
process: an anti-persistent FBM predicts that a merging 
step is anti-correlated with its immediate early occurrence 
of merging, and persistent FBM models positive correlation. 



3 SPHERICAL COLLAPSE 
3.1 Conditional mass function 

It is quite straightforward to solve the diffusion equation 
in spherical collapse model , we can actually copy the solu- 
tion in the literatures fe.g. iBond et alj [l 991: La cev &i Cold 
1 19931 ). Spherical collapse is equivalent to employ the bound- 
ary condition that there is a fixed absorbing barrier of height 
5 C = u)i — u>o > to the reformed random walk described 
by Eq. (0. The number of trajectories of the reformed walk 
within (S, 5 + dS) at S is 
1 



Q a d5 



-(8-26 c y/2S 



dS 



(10) 



and the number of trajectories absorbed by barrier within 
(S, S + dS) is given by 



f(S,5 c )dS = -dS-^L f° Q a d5 

dS i-oo 

5 S- 3/2 expf^ ] 
2tt \ 2S 



(11) 



which directly yields the universal conditional halo mass 
function 



/(Si 



- So, wi 
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■ exp 
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(12) 



v^FCSi-So)^ 1 "V 2(Si -Sp) 2 ", 

If a = 1/2 we recover the Eq. (2.15) in lLacev fc Cole] (|l993h . 

The conditional mass function is in fact more funda- 
mental than the mass function, since the mass function can 
be recovered from the co nditional m ass function by setting 
the limit So -> 0,w -> (|Panll2007l '). 

2a S c ( 5 2 r 



f(S,S c )dS 



■ exp 



2S 2c 



dS 



(13) 



727S-+ 1 

Thus, prompted by the relation between mass function 
and conditional mass function, we recognize that there is a 
systematical effect due to the finite volume of simulations. 
This can lead to an underestimate of the halo mass function 
in the large mass regime (more details are discussed in Ap- 
pendix A) . Recall that the halo model tends to over-estimate 
the amplitude of the three point correlation function of dark 
matter, which can be corrected (at least partly) by applying 
certa in arbitrary high mass cut-off to the halo mass func- 
tion (|Wang et all 12004 iFosalba et alj|2005h . To predict the 
three point correlation function of a simulation, a better ap- 
proach would be to use the halo mass function of Eq. ()A1|| 
to include the finite volume effects. 



3.2 Merger rate 

Now we will conside r the merger rate function in line with 
lLacev fc Cole] i| 19931 ). For a trajectory which has experienced 
a first up-crossing over uji at Si , the conditional probability 
of having a first up-crossing over u}2{^2 < wi) at S2(S2 < Si) 
in the interval dS2 is given by the Bayes formula, 

/(Sl,Ul|S2, W 2 )dSl/(S2, LU 2 )dS2 



/(S2, W2|Sl, U)l)dS2 



2a LU2(ui — (^2) 



12-k 
x exp 
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f(Si,u>i)dSi 
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2 \{Si- 

(ui > U2 , Si > S2) • 



S 2 ) 



- V2)~ , ^2 

S 2 ) 2q S| q 



S 2 (Si 

2 



si a 



(14) 

This function is interpreted as the probability that a halo of 
mass Mi at time t will merge to build a halo of mass between 
M2 and M2 + dM.2 at time t2 > ti. The mean transition 
rate is obtained from it by setting £2 — > ti (equivalently 

UJ2 — > OJl — Ul), 

1 a+1 



d 2 P {Si^s 2 n 



dS2dui 



rfS 2 dw = 
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Figure 2. Halo merger rate (left panel) and accretion rate (right panel) computed with Eq. 11161 . The accretion rate is simply AM /Mi X 
d 2 p/d\nAM/d\nt. Hereafter M* = 1O 13 M . 



Therefore the merger rate, i.e. the probability that a halo of 
mass Mi accretes a clump of mass AM — M 2 — Mi within 
time dixit (corresponding to dui), is 



than Zf, can be acquired by 



d 2 p(Mi -> M 2 \t)/dlnAMdhxt 



2a(M 2 )AM 



duo 
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dw 
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(16) 



In Figure [2] the halo merger rates of different masses 
as predicted by Eq. (|16p are plotted for comparison. Three 
models are presented: an anti-persistent FBM of a = 0.4, 
the standard normal Brownian motion of a = 0.5 and a 
persistent FBM of a = 0.6. 

For halo progenitors of the typical mass scale Mi = 
M*, the merger rate and halo accretion rate increases with 
increasing a at small AM /Mi, but at very large mass ratios 
(AM /Mi > 25) the case is reversed. 

This is also true for more massive progenitor masses 
(Mi = 10M*), but here the transition point in the mass 
ratio is smaller: AM /Mi = 3. 

For the less massive progenitors, Mi = 0.1M* , the tran- 
sition point is very high, that in the whole plotted range (up 
to AM /Mi — 100), the models with greater Hurst exponent 
always have greater merger and accretion rate. 



3.3 Halo formation time distribution 

The formation time (redshift) of a halo is the time (redshift) 
when its major progenitor contains half of the halo mass. 
According to the counter argument of LC93, the cumulative 
distribution of halo formation time, i.e. the probability that 
a halo of mass Mo at redshift zq is formed at redshift larger 



r s h M n 

P(>zf) = j s jjJ^y/(5l,*c(Z/)|5o,«e(«)))dSl (17) 



where So = S(M ) and S h = S{M /2). The halo forma- 
tion time distribution is simply given by the differentiation 
-dP(z f )/dz f . 

Figure [3] shows the predictions of the halo forma- 
tion redshift distribution for difference masses. In general 
smaller halos have more extended formation time distribu- 
ti on than th a t of la rger halos, in agreement with the results 
of iLin et ail (|2003l ). For a > 1/2, the halo formation red- 
shift distribution is more concentrated than the a — 1/2 
case, and for a < 1/2 the halo formation redshift distribu- 
tion is more extended. This seems to be in accordance with 
our finding that the merger rate is greater for larger a. 

As can be seen, the detailed effects of incorporating cor- 
relation between merging steps are intricate, varying with 
the halo mass. Impact on the halo formation time distribu- 
tion is small for those halos with masses around M*, and 
it becomes apparent only when halo mass deviates signifi- 
cantly from M*. The impact is also different for halos with 
small masses and large masses. Compared with the results 
given by the standard excursion set theory (a = 1/2), 

(i) M M* : negative correlation (a < 1/2) shifts the 
halo formation time distribution curve to the side of smaller 
Zf, which means that the younger halos are more abundant 
and the older ones are less abundant; whilst for positive 
correlation (a > 1/2) the older ones are more abundant; 

(ii) M > M* : the impact of the correlation is opposite to 
the case of small halo mass, negative correlation (a < 1/2) 
boosts more halos to form at an earlier time while positive 
correlation (a > 1 /2) induces more halos to form at an later 
time. 
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Figure 3. Halo formation time distributions for halos at z = with mass of 0.1, 1 and 10M*, which are predicted by the modified 
excursion set theory of different a assuming spherical collapse. The left panel is the cumulative distribution. 



4 ELLIPSOIDAL COLLAPSE 



and 



4.1 Moving barriers 

The spherical collapse model is perhaps accurate at high 
redshift e.g. the re-ionization era. However, as shown by sim- 
ulations, at low redshift the collapse of a clump of mass is 
ellipsoidal. For the excursion set theory, the general collapse 
condition for halo formation is not a constant 8 C any more, 
but a moving barrier B(S). 

Imposing a moving barrier on the random walk of 
(S, AS) described by Eq. © is equivalent to the case of 
placing a constant barrier boundary cond ition to the d iffu- 
sion equation with an extra drifting term (|ZentnerJl2007l ) . In 
our case of FBM, the Fokker-Planck equation turns out to 
be 



dQ a 
OAS 



= aAS z 



i d 2 Q a , dAB dQ 



OS 2 



+ 



dAS 05 



(18) 



where AB = B(Si) — B(So)- Unfortunately this equation can 
be solved analytically only for a few very special cases (see 
Appendix B). 

4.2 Conditional mass function: the Sheth-Tormen 
approximation 

Although we do not have analytical solution of Eq. (|18l) . it 
could be solved with the numerical method of lZhang fc Huil 
(|2006l ). However, this is inconvenient for general ex- 
plorations. Here we b orrow the pragmatic approach of 
ISheth fc Tormenl (|2002l ). the conditional mass function is 
then approximated by 

2a \T(Si, zi\S , zo)\ 



/(Si - S )dSi = dSi- 



x exp 



/2¥ (Si 
[B(Si,zi) 



So) Q+1 
-B(S ,z )} 2 



2(51 - So) 2 



in which 



T(Si|S ) = £ 



(So - Si)" n [B(S!) - B(S )] 
n\ dSi n 



(19) 



(20) 



B(S,z) = ^S c (z) 



1 + 6 



qSi(z) 



(21) 



with q = 0.707, b = 0.4 85 and 7 = 0.615 (|Sheth et alj|200ll : 
ISheth fc Tormenll2002l ). 

Inspired by the results given in Appendix B, it is prob- 
ably more appropriate to use 



(-S)" d n AB 
n! dS n 



(22) 



with S = (Si — So) 2a in Eq. 1191 Testing the two approxima- 
tion methods numerically, we found that the actual differ- 
ence between the two is in fact very small. So we opted to 
use the numerically simpler Eq. (|20[) . 

In figure|?]we compare our prediction on the conditional 
mass fu nction with that m easured from the Millennium Sim- 
ulation (|Cole et alJboOSl ) at z = 0,zi = 0.5,1,2,4, for 
halos of masses Mo=0.1, 3.16, 100A/*. It is indeed very 
interesting to notice that the performance of the origi- 
nal Sheth-Tormen approach (Eq. [19] of a = 1/2) is fairly 
good though no t as remarkable as t he fitting functio n of 
ICole et al.l (|200cf ). It is surprising that lCole et al.l (120081 ) ex- 
trapolated the Sheth-Tormen mass function for the condi- 
tional mass function in ellips odal collapse model, in s pite of 
the one actually proposed in lSheth fc Tormenl (|2002T l . 

From Figure [4] we can see that the effects of correla- 
tion between merging steps gradually decrease with time, 
with halo formation becoming sensitive to the correlation at 
high redshift. Except in cases of high redshifts and of ex- 
tremely high halo mass, it appears that positive correlation 
(a > 1/2) is preferred by the simulation in ACDM universe. 
Nevertheless, one must be cautious to this result, as the mea- 
surement presented is rather crude and lacks error bars. Its 
accuracy is not sufficient to justify with confidence whether 



2 The data is publicly available 

http: / '/star-www. dur. ac.uk/~cole/merger_trees 



at 
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Figure 4. Conditional mass functions of progenitor halos at rcdshift z\ =0.5, 1, 2, 4 (identical for each row) for t hree different hal o 
masses Mo (each column has the same mass). Histograms are the results from the Millennium Simulation given bv lCole et al ] l|2008l) . 
while lines are predictions by Eq. H19H of three different values of a as indicated. 



the correlation between merging steps is positive, negative 
or zero. 

4.3 Halo formation time distribution 

The halo formation time distribution is recalculated with the 
modified excursion set theory in ellipsoidal collapse model 
and the results displayed in Figure [5] Although the halo for- 
mation time distributions are modulated significantly by the 
non-constant barrier, its dependence on a does not change, 
and demonstrates similar trends as the spherical models pre- 
dicts: positive correlation results in more old small mass ha- 
los and more young large mass halos, while the effects of the 
neg ative correlation are to the contrary. 

iGiocoli et af] (|2007l ) calculated the cumulative distri- 
bution of halo formations time and compared with GIF2 
simulation. They found that the median formation redshift 
of halos of mass less than M* extracted from the simula- 
tion is larger than that predicted by the Sheth-Tormen for- 
mula, even after some numerical improvement. Combined 



with our discovery of effects of correlation, it seems that the 
discrepancy can be explained by the existence of a positive 
correlation between adjacent merging steps. 



5 CONCLUSIONS AND DISCUSSION 

In this work the FBM is formally incorporated in the excur- 
sion set theory, to account for possible correlations between 
merging steps of halo formation. Such modification is min- 
imal but provides a unified theoretical frame to investigate 
the effects of correlation. More specifically, the correlation 
depends on the Hurst exponent a, for a < 1/2, the correla- 
tion between adjacent merging steps is negative, for a > 1/2 
the correlation is positive, while a = 1/2 is reduced to the 
standard case of null correlation. 

We calculated the conditional mass functions in models 
of the spherical collapse and the ellipsoidal collapse. Differ- 
ent collapsing conditions do not change the relative effects 
of correlations. It is revealed that for a typical progenitor of 
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Figure 5. Halo formation time distributions for halos of masses 0.1, 1 and 10M* at z = 0, which are predicted by the modified excursion 
set theory of different a in ellipsoidal collapse model. The left panel is the cumulative distribution. 



M* , a positive correlation will boost the merger rate when 
the mass to accrete is less than ~ 25M*, and negative cor- 
relation will significantly reduce such a rate. Accordingly we 
find that compared with the standard excursion set theory 
without correlation, positive correlation increases the aged 
population for small mass halos but produces more young 
members for large mass halos; while the negative correlation 
projects opposite effects to the halo formation. The same 
trend is also seen for other progenitor masses, although the 
transition mass ratio is greater for small progenitor masses 
and smaller for greater progenitor masses. 

The results of the modified excursion set theory in ellip- 
soidal collapse model are checked w i th tho se me asurements 
from simulations of iGiocoli et al.l |2007l ) and ICole et all 
(|2008l ) . The comparison indicates that there is a sign of weak 
positive correlation, although not yet of significant confi- 
dence, as it is limited by the accuracy and dynamic ranges 
of their analysis. The c onclu sion appears to be in conflict 
with the claim of IPanl (|2007l '). but note that in that work 
the mass function was calculated with the spherical collapse 
model, which is not valid a t mo s t reds hifts of concern. Be- 
sides, the main purpose of IPanl (|2007l ) is to show that the 
effect of correlation is not trivial. 

The accuracy of the Sheth-Tormen approximation to 
the ellipsoidal collapse in the modified (FBM) excursion set 
theory has not been checked against rigorous numerical so- 
lutions. However, since the actual correlation is weak, i.e. 
a only deviates from 1/2 slightly, the approximation suf- 
fices for general qualitative discussion. Calibration and sub- 
sequent improvement of the approximation with numerical 
computation is of course needed for more precise analytical 
modeling and implementation in semi-analytical models of 
galaxy formation. 

In this paper we explored the effects of correlations be- 
tween mergers on the global halo formation process. It is 
of interests to actually check if such a correlat i on wil l lead 
to the assembly bias discovered by iGao et all (120051 ) , and 
whether it can explain the dependence of halo formation on 
the large scale environment. 



A by product of our work is the discovery of a system- 
atic effect in the mass function due to the finite volume of 
simulation. It explains why applying a high mass cut-off to 
the mass function can improve the halo model prediction on 
three point correlation function. A better approach in the 
halo model calculation of three point function is to shift the 
theoretical mass function according to the formulae given 
in Appendix A, this may lead to better agreement between 
simulation and analytical results. 
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APPENDIX A: FINITE VOLUME EFFECTS IN 
SMALL-BOX SIMULATIONS 

N-body simulations are always performed in cubic boxes of 
finite size L s ; m . By definition, the density contrast of whole 
the simulation box is zero, but the variance is not. Hence- 
forth the measured mass function of simulation is condi- 
tional, determined by /(Si = S, Wi = 5 c \So{L sim ) ^ 0, u)q = 
0). From Eq. ()12p the explicit halo mass function of a sim- 
ulation is then 



f(<7, L aim )dh\a 



4a 



where 



and 



x exp 



S = a = (2tt) 



2tt (° 2 - 



2(o 2 -* 2 ) 



a 2 ) a a 2 - ( 
din a 
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So = o~ (L s 
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f oo 

7 
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The infrared cutoff in the integration for a 2 accounts for the 
deficiency of large scale power in finite box, w hich has al- 
ready been identified as a systematic effect (e.g. iReed et al.l 
l2007f ). 

The difference between Eq. (|A1[I and Eq. (|13[) is shown 
in Fig. IA1I for different box sizes. The fiducial model has 
a = 0.436, by which the spherical model provides good ap- 
proximation to_the mass function of simulations in high mass 
region i|Panll2007T ). We can see that the finite volume effect 
becomes apparent only when — In a > 0.2, i.e. in high mass 
regime, and the resulting fractional error quickly drops down 
to < 10% for simulations of L 3 i m > 100/i _1 Mpc. 

Nowadays it is a very common practice to use simu- 
lations with boxes as small as < 50/i _1 Mpc to extend the 
dynamic range. The finite volume effect demonstrated in 
Fig. IA1I has to be taken into account for precision mea- 
surement. The necessary correction is actually very simple: 
the true measured mass function f(S)dS can be recovered 
by shifting the raw mass function f(S\So)dS along S by 
So. Note that there is a cosmic variance for So, for each 



Figure Al. Deviation of mass functions in finite-sized simula- 
tions (Eq. lAH to Eq. lT3t . For the fiducial model a = 0.436. -L s ; m 
is the size of the cubic box for simulation, ao is calculated with 
infrared cutoff in integration (see text), and a spherical top- hat 
window function in real space of radius L s ; m /(47r/3) 1//3 . Devia- 
tion is negligible in regime — lncr < 0.2. 



individual realization of simulation one m ay need to calcu- 
late the individua l So from the realization (|Reed et alj|2007l ; 
iLukic et al1l2007l ). 



APPENDIX B: PSEUDO-LINEAR BARRIER: 

B = 5 C + f3S 2a 

Starting with Eq. <jTHJ) , let So = and Si = S > 0, a variable 
substitution S = S 2a yields 

dQ a 1 d 2 Q a 



dS 



1 dB dQ a 
2 dS 2 2aS 2a ~ 1 dS Q5 



If the barrier is B = 5 C + f3S 2a , the above nonlinear Fokker- 
Planck equation is further simplified to 

dQa 1 d 2 C 



dS 



2 OS 2 85 



(B2) 



which is the well known pr oblem of a linear barrier up on the 
normal Brownian motion (|Shethlll998l ; IZentnenl2007T l . 

The probability of a trajectory that has the first barrier 
passage within (S, S + dS) in the solution is given by 



f(S)dS = 



= — exp 



2ttS 3 / 2 



{PS + 5_c 
2S 



dS , 
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from which the universal mass function is easily obtained 



/((j)dlncr = 



4a 5 C 



2tt cr 2 



■ exp 



(5c + l3a A 
2ct 4q 
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However, there is no analogous analytical expression for 
the conditional mass function under this type of barrier, as 
the barrier does not have the symmetry that a true linear 
barrier possesses: for any pair of points at So and Si > So 



B 1 -Bo=P{Sl a -S^)^^ 



So) 



(B5) 



This is why we call it pseudo-linear. 



